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Abstract 

Using a large sample (~ 11800 events) of — > i^^vr^e+z/g and —t- K~tt~^ fi'^u^^ decays 

collected by the CLEO-c detector running at the ^^(3770), we measure the helicity basis form 

factors free from the assumptions of spectroscopic pole dominance and provide new, accurate 

measurements of the absolute branching fractions for K e^Ug and — > K fi'^u^ decays. 

We find branching fractions which arc consistent with previous world averages. Our measured 

helicity basis form factors are consistent with the spectroscopic pole dominance predictions for 

the three main helicity basis form factors describing K*^l'^ui decay. The ability to analyze 

-D+ K~Tr~^ I^^Vfj^ allows us to make the first non-parametric measurements of the mass-suppressed 

form factor. Our result is inconsistent with existing Lattice QCD calculations. Finally, we measure 

the form factor that controls non-rcsonant s-wave interference with the K*^l'^ve amplitude 

— *o 

and search for evidence of possible additional non-resonant d- or /-wave interference with the K . 
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I. INTRODUCTION 



We present new measurements of the — )■ K e^v^ and — )■ K yU^z/^ absolute 
branching fractions, their ratio, and measurements of the semileptonic form factors control- 
hng these decays.^ Exclusive charm semileptonic decays provide particularly simple tests 
of over decay dynamics since long distance effects only enter through the hadronic form 
factors [l| . A wide variety of theoretical methods have been brought to bear on the calcula- 
tion of these form factors including quark models [2j], QCD sum rules Lattice QCD 
analyticity jsf, and others j^. Using a technique developed by FOCUS [3], we present non- 
parametric measurements of the dependence of the helicity basis form factors that give an 
amplitude for the K't:^ system to be in any one of its possible angular momentum states 
where is the invariant mass squared of the lepton pair in the decay. The ultimate goal of 
this study is to obtain a better understanding of the semileptonic decay intensity. 

CLEO-c produces D mesons at the '?/'(3770), which ensures a pure DD final state with no 
additional final state hadrons. In events where the — )■ K'n^i^ui is produced against a 
fully reconstructed D~ the missing neutrino can be reconstructed with unparalleled precision 
using energy-momentum balance. Hence, CLEO-c data offer unparalleled q^ and decay angle 
resolution allowing one to resolve fine details in the structure of these form factors without 
the complications of a deconvolution procedure. The various helicity basis form factors are 
distinguished based on their contributions to the decay angular distribution. 
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FIG. 1: Definition of the Oy, Og, and x angles. 



The amplitude A for the semileptonic decay — )■ K'n'^i'^i'i is described by five kine- 
matic quantities: g^; the kaon-pion mass {rriKn)] the kaon helicity angle (^y), which is 
computed as the angle between the vr and the D direction in the K~7i^ rest frame; the lep- 
ton helicity angle {6e), which is computed as the angle between the ui and the D direction in 
the i'^h'e rest frame; and the acoplanarity angle between the two decay planes (x)- The decay 
angles are illustrated in Fig. [H The amplitude A can be expressed in terms of four helicity 
amplitudes representing the transition to the vector K*^: H^lq"^), H_{q^), Hoiq"^), Ht{q^) 



^ Throughout this paper the charge conjugate is imphed when a decay mode of a specific charge is stated. 
^ We reconstruct _D+ K*^ t^v^ modes as K^n^ l^v^ decays, and use the Clebsch-Gordan factor 

1.5 to correct for K K decays, which we do not detect. 
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and a fifth form factor, ho{q'^)HQ{q'^) describing a non-resonant, s-wave — t- K~7T~^i~^i>i 
contribution. 

The diferential decay width for the 4-body semileptonic process is 



d cos 9i dcosOy dx dq^ dm\^^ 25QTC^m'l)y/q^mK-n 



(1) 



where is the decay intensity, K is the K tt^ momentum in the rest frame, P* is 
the momentum of the kaon in the K^tt^ rest frame, and |P^| is the momentum of the in 
the i^u rest frame. Upon integration over x, the differential decay width is proportional to: 



((l + cos^,) smey)'\H4q')m' 
/l,|2, _ q'-mj ) +((l-cos^,)sinev)'l^-(5')l'l/^l' 

+8 sin^ eecoseyHo{q'^)ho{q'^)Re{Ae-'^ (3} 



(2) 



(sin sin ^V)' I ^+(9^)1^ + {smeesineyyiH^q^^^^ 



+ ^{q' - miY^ I +{2 cos 9^ cos ev)'l^o(g') P 

1 [ +{2coseyf\Ht{q')\^ + 8coseeCOS^eyHo{q^)Ht{q^) 



The Ht{q^) form factor, which appears in the second term of Eq. (|2]), is hehcity suppressed 
by a factor of vn\jq^. The mass-suppressed terms are neghgible for — K'-n^e^Ve but 
can be measured in — )■ i^^vr+yU+z/^. The Ht{q^) form factor can only be effectively 
measured in — > K^-k^ ^^v^ decays at low q^ where the mass suppression effects are least 
severe. The semimuonic to semielectric branching ratio is sensitive to the magnitude of the 
Ht{q^) form factor. 

We study the form factor of the non-resonant, spin zero, s-wave component to — >■ 
K*^ li^v^ first described in Ref. [S'j. According to the model of Ref. joj, 2.4% of the decays 



in the mass range 0.8 GeV/c < rriK-K < 1-0 GeV/c are due to this s-wave component |15 



where mKn is the K tt~^ mass. The underlined term in Eq. ([2]) represents the interference 
between the s-wave, K~7r~^ ampliti 
Breit-Wigner function of the form: 



between the s-wave, K~Tr~^ amplitude and the K amplitude, represented as a simplified. 



/3 = ^ — (3) 

mj^^ -ml + imoT (^) 

where P* is the kaon momentum in the K^tt^ rest frame, and Pq is the value of P* when 
the K~7r~^ mass is equal to the K* mass^. 

The s-wave form factor is denoted as holq"^) in the underlined piece of Eq. ([2]). Following 
Ref. [3] we model the s-wave contribution as an amplitude with a phase {6) and modulus 



^ We are using a p-wave Breit-Wigner form with a width proportional to the cube of the kaon momentum 
in the kaon-pion rest frame. Our Breit-Wigner intensity is proportional to P*^ as expected for a p-wave 
Breit-Wigner resonance. Two powers of P* come explicitly from the P* in the numerator of the amplitude 
and one power arises from the 4-body phase space as shown in Eq. ([1]). We are not including additional, 
small corrections such as the Blatt-Weisskopf barrier penetration factor. 
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(A) that are independent of mK-K- We have dropped the second-order, s-wave intensity 
contribution (oc in Eq. (|2]) since A <^ \/3\. 

The X integration significantly simphfies the intensity by ehminating all interference terms 
between different helicity states of the virtual with relatively little loss in form factor 
information. 

The four helicity basis form factors for the — )■ K /i"*"//^ component are generally 
written fToJ as linear combinations of a vector {Vlq"^)) and three axial- vector (Ai 2,3(0'^)) 
form factors according to 



[Ml - - g2)(Mz5 + mK^)A^{q^) 



2mKnVq^ 



MdK 



-4- 



(4) 



{MD + mK^)Ai{q^) - 



{Ml - + g') 
Md + rriK-K 



2g2 



Md + rriKn 



Asiq') 



where Md is the mass of the and K is the momentum of the K~7r ' sy 
frame of the D~^. In the Spectroscopic Pole Dominance (SPD) model (ol. llo|. 
vector form factors are given by 



system in the rest 
these axial and 



Viq 



V{0) 



1 



A,2,3(0) 

- - qVMl 



(5) 



where My = 2.1 GeV/c^ and Ma = 2.5 GeV/c^. The SPD model allows one to parameterize 
the H_{q^), H^{q^), HQ^q'^), and Ht{q^) form factors using just three parameters, which 
are ratios of form factors taken at = : ry = ^(O)Mi(O), r2 = A2{0)/Ai{0) and 
r3 = A3{0)/Ai{0). There are accurate measurements [91] of ry and r2, but very little is 
known about r^, which is an important motivation for this work. 

In this paper, we use a projective weighting technique [3] to disentangle and directly mea- 
sure the dependence of these helicity basis form factors free from parameterization. We 
provide information on the six form factor products if^(g^), ifQ(g^), ho{q'^) HQ{q'^) , H^{q^) 
and HQ{q^)Ht{q^) in bins of q^ by projecting out the associated angular factors given by 
Eq. ([2]). We next describe some of the experimental and analysis details used for these 
measurements. 



II. EXPERIMENTAL AND ANALYSIS DETAILS 



The CLEO-c detector [llj consists of a six-layer inner stereo-wire drift chamber, a 47- 
layer central drift chamber, a ring- imaging Cerenkov detector (RICH), and a cesium iodide 
electromagnetic calorimeter inside a superconducting solenoidal magnet providing a 1.0 T 
magnetic field. The tracking chambers and the electromagnetic calorimeter cover 93% of the 
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full solid angle. The solid angle coverage for the RICH detector is 80% of An. Identification 
of the charged pions and kaons is based on measurements of specific ionization (dE/dx) in 
the main drift chamber and RICH information. Electrons are identified using the ratio of 
the energy deposited in the electromagnetic calorimeter to the measured track momentum 
{E/p) as well as dE/dx and RICH information. Although there is a muon detector in CLEO, 
it was optimized for b-meson semileptonic decay, and is ineffectual for charm semileptonic 
decay since a muon from charm particle decay will typically range out in the first layer of 
iron in the muon shield. 

In this paper, we use 818 pb~^ of data taken at the V'(3770) center-of-mass energy with 
the CLEO-c detector at the Cornell Electron Storage Ring (CESR) e^e~ collider, which 
corresponds to a (produced) sample of 1.8 million D^D~ pair events [l^ . 

We select the events containing a D~ decaying into one of the following six decay modes: 
D- K'^s^-, D' K+n-n^, K°sn-n°, K+n-n-n°, K^gT^- -r- 7:+ , and 

— 7- K~ K^T{~ along with a 4-body semileptonic candidate. To avoid complications due 
to having two or more — )■ K^n^i'^i'i decay candidates in the event, we select the decay 
candidate with the smallest \M-^^^ — Mj^- \ value where M-^^ is the beam-constrained mass. 



The beam-constrained mass Mj^^ is defined as M-^^c? = y K-E'beam)^ ~ c^-^qI where E-^^^^^ 
is the beam energy and Pd is the D-tag momentum. More details on selecting the tagging 



candidates as well as identifying 7r° and Kg candidates are described in Ref. 12 . 

We used extensive Monte Carlo (MC) studies to design efficient, background-suppressing 
selections. The — )■ K'n'^i'^ui reconstruction starts by requiring three well-measured 
tracks not associated with the tagging D~ decay. In order to select semileptonic decays, we 
require a minimal missing momentum and energy of 50 MeV/c and 50 MeV, respectively. 
Both the minimal missing momentum and energy are calculated using the center-of-mass 
momentum and energy. In order to reduce backgrounds from charm decays with missing 
7r° 's, we require an unassociated shower energy of less than 250 MeV. The unassociated 
shower energy refers to electromagnetic showers, which are statistically separated from all 
measured, charged tracks. Charged kaons and pions are required to have momenta of at 
least 50 MeV/c and are identified using dE/dx and RICH information. We require that the 
pion deposits a shower energy, which is inconsistent with the electron hypothesis. 

Electron candidates are required to have momenta of at least 200 MeV/c, lie in the 
good shower containment region (|cos^| < 0.9), and pass a requirement on a likelihood 
variable that combines E/p, dE/dx, and RICH information. Our simulations indicate that 
contamination of our kaon sample due to pions is less than 0.06% using this likelihood 
variable. The only final state particle not detected is the neutrino in the semileptonic decay. 
The neutrino four-momentum vector can be reconstructed from the missing energy and 
momentum in the event. The resolution, predicted by our Monte Carlo simulation, is 
roughly Gaussian with an r.m.s. width of 0.02 GeV^/c^, which is negligible on the scale that 
we will bin our data. 

For — )■ K~'ir~^ fi^Ufj^ candidates, it is difficult to distinguish the tt^ track from the 
track. Because — )■ f^^vr+zi+z/^ decay is strongly dominated by — )• K~n~^, which is 
a relatively narrow resonance, we select the positive track with the smallest ImKn — fTT'j^*o\ 
as the pion and the other track as the muon. Our Monte Carlo studies concluded that this 
K arbitration approach was correct 84% of the time and works better than pion-muon 
discrimination based on the electromagnetic calorimeter response. 

We apply a variety of additional requirements to suppress backgrounds in — )■ 
K^Ti^ fi^Un candidates. We require that the muon is inconsistent with the electron hy- 
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pothesis according to the electron likelihood variable. We require that missing momentum 
(-^miss) ^^^^ within 20 MeV of the missing energy (-Ej^jgg). For _D+ — )• K'n^ [I'^i/^j, candi- 
dates, we also require -0.01 < M2 < 0.015 GeVVc^ The distributions for muon 
and -0.01 < < 0.015 GeVVc^ The distributions for muon and electron 

I 1 lliloo 

candidates are illustrated in Fig. [21 

In order to suppress cross-feed from — )■ K'-K'^e^Ve decay to our — ?■ K^-k^ ^'^v^ 
sample, we construct the squared invariant mass of the lepton candidate, M^c^ = 

{^E-^Q^^ — -E]3^ag — — Ek — -Ett) — (cP^)^, where E-q^^^ is the reconstructed energy 
of the D~ produced against the — )• i^^Tr+f+z/^ candidate and Ek, E^,, Pe are^ the 
reconstructed kaon energy, pion energy, and lepton momentum. We require < < 
0.020 GeVVc^ to eliminate both K'Ti+e+Ue cross-feed and D+ K-Tt+Tx+Ti^ de- 

cays. In order to suppress backgrounds to — )■ K^n^ fi'^u^^ from — )■ K^n^n^ decays 
with an accompanying bremsstrahlung photon, we require that cosine of the minimum angle 
between three charged tracks and the missing momentum direction be less than 0.90. This 
requirement is illustrated in Fig. |3l 

We obtain 11801 — i- K^ir^i^Ui candidates. The mKn distribution for these 
—7- K^TT^i^ue candidates is shown in Fig. |H Finally, we require 0.8 GeV/c^ < mKir < 
1.0 GeV/c^ and select 10865 events. 

Two types of Monte Carlo simulations are used throughout this analysis. The generic 
Monte Carlo simulation is a large charm Monte Carlo sample consisting of generic DD 
decays, which is primarily used in this analysis to simulate the properties of backgrounds 
to our — )■ K^TT^i^Ui signal states. The generic Monte Carlo events are generated by 
EvtGen 13] and the detector is simulated using a GEANT-based [13] program. In much 
of the form-factor work, we use an SPD Monte Carlo simulation based on the SPD model 
described in Sec. [land summarized by Eqs. ([2H5]). We use the SPD parameters of Ref. jof, 
ry =1.504 , r2 =0.875, and we set r3=0. 

The background shapes in Fig. [H are obtained using generic Monte Carlo simulations. 
Our simulation predicts a 6.5 % background for our — )■ K~tt^ ^x^v^ sample with 4% 
due to misidentified — > K^Ti^e'^Vf. cross-feed events and the rest due to various charm 
decays. The simulation also predicts a 1% background to our — )• K^n^e^Ve sample with 
0.03 % due to — )■ K^n^ fi'^u^ cross-feed. 



III. ABSOLUTE AND RELATIVE BRANCHING FRACTIONS 

We have measured both the semimuonic to semielectric relative branching ratio and the 
— )• K e'^Ve and — )■ K /i"*"//^ absolute branching fractions, which we will denote as 
and respectively. The B^/Be relative branching ratio is expected to be less than 1 due to 
the reduced phase space available to the semimuonic decay relative to the semielectric decay. 
The mass-suppressed terms in Eq. ([2]) will change the relative branching ratio compared to 
the phase space ratio. In the context of the SPD model, Eq. ([S]), the relative branching 
fraction will depend on = A3{0)/Ai{0), which controls the strength of the Ht{q^) form 
factor and is essentially unknown. It is expected that Bf^/Be will increase with increasing 
values of r^. 

In order to obtain the semimuonic to semielectric branching ratio, we write the observed 
D+ -> K*^fi+u^ and D+ K*°e+Ue yields as 
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M2. (GeV^/c") 

miss \ ' 

FIG. 2: The -^^jgg distributions for events satisfying our nominal — )■ K^ir^l^vg selection 
requirements apart from the ^^^gg requirement, (a) shows the ^^^gg distribution for — >■ 
K^TT^ ji^v^ candidates, while (b) shows the ^^jgg distribution for — )■ K^ir^e'^i'e candidates. 
For —7- K^t:'^ fi'^u^ candidates, we require that -^^jgg hes between the vertical lines. This cut 
is placed asymmetrically on our semimuonic sample to suppress cross-feed from — )• K~TT~^e~^Ue. 
In each plot, the solid histogram shows the signal plus background distribution predicted by our 
Monte Carlo simulation, while the dashed histogram shows the predicted background component. 





(6) 

where |/^,e are the observed yields, b^^e are non-semileptonic backgrounds, and n^^e give 
the number of produced semileptonic decays in our data sample. The cross-feed matrix, 
which multiplies the Ue and signal vector, is constructed from e^^eif)) which are the 
—7- K*^^^!/^ and — t- K*^e~^i'f. detection efficiencies, and c^,e(/); which are the cross- 
feed efficiencies. For example, c^(/) is the efficiency for reconstructing a — )■ ^*°/i+z/^ 
event as a — )■ K*'^e~^i'e candidate. The y^^e yields are obtained by counting the number 
of semimuonic and semielectric events in our mass range 0.8 < m^^^ < 1.0 GeV/c^. The 
relative branching ratio is given by B^/Be = n^/rie. 

The vector / represents parameters that the efficiencies and cross-feeds can depend on 
such as the SPD parameters: ry , ^"2, and and the s-wave amplitude and phase. The 
detection efficiencies, e^,e(/), and the cross-feed efficiencies, c^,e(/), were obtained using 
our Monte Carlo simulations. We will refer to the use of Eq. ([6]) to obtain the relative 
branching ratio, B^/Be-, as the cross-feed method. 



We used the double-tag technique, described in Ref. [12||, to measure the K*^jj^ 
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COS e 



FIG. 3: Distributions of the largest cosine between missing momentum vector and any of the three 
charged tracks from the semileptonic candidate when ah cuts are apphed except the cut on largest 
cosine, (a) shows the cos 9^^^ distribution for D'^ — t- K~ti~^^'^v^ candidates, while (b) shows 
the cos^™^^ distribution for — )• K~'K~^e'^Ve candidates. We remove all combinations to the 
right of the vertical line, which removes the major part of remaining Kttit background for the 
semimuonic sample. In each plot, the solid histogram shows the signal plus background distribu- 
tion predicted by our Monte Carlo simulation, while the dashed histogram shows the predicted 
background component. 

and — )■ K*^e^Vf. absolute semileptonic branching fractions (-B^,e)- We define single tag 
(ST) events as events where the D~ was fully reconstructed against one of our six tag modes 
without any requirement on the recoil . 

We estimate the number of ST events by fitting the distributions, shown in Fig. [5l 
using a binned maximum likelihood fif^. Here AE = Ej^ — E^^^^, where Ej^ is the energy 
of the D-tag candidate. 

The total number of reconstructed ST events is then 

"-ST ^ ^^+£'-egxStag' (''') 

where 'T-g'j is the number of ST reconstructed events in the i-th mode, ND+ is the number 
of produced D^D^ events in our data sample, egrj is the ST detection efficiency, and -B^g^g 
is the tag mode branching fraction. 

For double tag (DT) events, we reconstruct into one of our six tagging modes, and 
require the presence of either a — )■ K*^ ^'^v^ or — t- K*^e^Ve candidate. The DT yields 



^ Our fitting function is a sum of Gaussian and Crystal Ball line-shape functions [12] over a first order 
Chebyshev background polynomial. 
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FIG. 4: The mx-K distributions for events satisfying our nominal — t- K^ir^i^i'i selection re- 
quirements, (a) shows the itik-k distribution for — t- K^tt^ ^'^u^ candidates, while (b) shows the 
rriK-K distribution for — t- K~Ti~^e^Ve candidates. Over the full displayed mass range, there are 
11801 (6 227 semielectric and 5 574 semimuonic) events satisfying our nominal selection. For this 
analysis, we use a restricted mass range from 0.8 - 1.0 GeV/c^ which is the region between the ver- 
tical lines. In each plot, the solid histogram shows the signal plus background distribution predicted 
by our Monte Carlo simulation, while the dashed histogram shows the predicted background com- 
ponent. In this restricted region, there are 10865 ( 5 658 semielectric and 5 207 semimuonic) events. 
The inserted figures are on a finer scale to better show the estimated background contributions. 



are then 

= Nj,,n- (e^T^tag^e + cf^V^tag^.) (») 

and 

^DT = ^D^D" (^DT^tag^e + e^'r^iS^ag^^) ' (9) 

respectively. The yields ^j^*rj and ^£)*rj represent the number of reconstructed DT events in 
semielectric and semimuonic decay modes after the background subtraction. The efficiencies 
eg*rj and event detection efficiencies for the semielectric and semimuonic 

decay modes. The cross-feed efficiency cj^rj describes how often a semimuonic decay is 

reconstructed as a semielectric candidate, while the cross- feed efficiency cf^j, describes how 
often a semielectric decay is reconstructed as semimuonic candidate. The variables B^, 
are the respective K*^e^Uf. and K*^^'^v^ branching fractions, which we wish 

to measure. 

Dividing Eq. (ED and Eq. ^ by Eq. (El), we have: 

(10) 
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Equation (fTOl) shows how the branching fractions of — )■ K*^e^Ve and — )■ K*^fi~^i'n 
semileptonic modes depend on the ratio of the DT and the ST yields, the detection efficien- 
cies, and the cross-feed efficiencies. 

Figures O and [7] shows the AE distributions for our double tag sample. For both semilep- 
tonic decay modes, about half of our sample comes from the — t- K^ti^h^ D-tag mode. 
The ST yields for this mode are nearly background free. The cross-feed fraction for the 
D~ — )• K*^e~i'e semileptonic mode is less than 0.02%, while, for the D~ — >■ K^'^n'u^ semilep- 
tonic mode, the cross-feed fraction is 3.7%. The background level is about 2.5 times smaller 
for the D'^ — )■ K*^e~Ue mode than for the — t- K*^ ^"v^ mode. The semielectric mode is 
nearly background free because of the effectiveness of the electromagnetic calorimeter, while 
our semimuonic mode uses a variety of less effective kinematic cuts to suppress background 
and cross-feed. 

Our absolute branching fraction results are summarized by Tables [I] and [ITl Table [T] gives 
a "conditional" absolute branching fraction based only on — ?■ K^tt'^I^ui decays into the 
mass range 0.8 < rnKv < 1-0 GeV/c^. This mass range is required for events entering into 
Figs. [HI and [71 We find that the total systematic error for the semielectric and semimuonic 
absolute branching fractions, presented in Table I, are comparable. The dominant system- 
atic error for the semielectric decay is due to the 1% uncertainty in the efficiency our electron 
identification requirements, while the dominant systematic error for the semimuonic branch- 
ing fraction is due to the 0.8% uncertainty in the background subtraction. The remaining 
systematic error, which is 1.2% for both the semielectric and semimuonic branching frac- 
tions, includes uncertainties in the final state radiation corrections, as well as uncertainties 

in the tracking and particle identification efficiencies for the kaon and pion tracks. Table [TTl 

— *o 

on the other hand, relies on models for the K line-shape to extrapolate outside of the 200 
MeV/c^ wide mass region where our measurements are made in order to report the conven- 

tional — )• K i'^u^ absolute branching ratios, which includes events over the entire rriKv 
spectrum. We include an additional, Clebsch-Gordan factor of 1.5 in order to correct for 
the undetected K* K tt^ decay mode^. Finally, we have included an additional ±0.10 % 
contribution to the quoted systematic error in Table [TTl based on the difference between the 
K extrapolations made using our Generic and SPD models. This ±0.10 % systematic 
error contribution includes both distortions to the K line shape as well as uncertainties in 
level of non-resonant contributions due to the s-wave amplitude. 

TABLE I: Conditional absolute branching fractions. These branching fractions only represent the 
K^TT^ spectrum from 0.8 < tjik-k < 1-0 GeV/c^. 





Mode 




Branching fraction [%] 


r{D+ ■ 




-^e)/r(Z)+) 


3.19 ±0.04 ±0.05 


r{D+ - 






3.05 ±0.04 ±0.05 



Figure [8] and Table [TTT] compare our relative B^/ obtained using the cross-feed method 
to the ratio of absolute branching ratios for the six tag states and generic and SPD Monte 



^ The central values reported in Table HIl assume that all of the signal events in the 0.8 < rriKTr < 1-0 GeV/c^ 
mass region, where our AE measurements made, are due to Z?^ — t- K I'^i'i decay. 
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TABLE II: Comparison of our absolute branching fraction measurements to previously published 

data. These branching fractions represent the Kt: contribution over the full rriK-K spectrum and 

— *o 

include a systematic error contribution for uncertainties in the K line shape. 







Lumin. [pb ^] Be [* 










These results 




818 


5.52 


± 


0.07 




0.13 


CLEO [15] 




56 


5.56 




0.27 




0.23 


World Average 


[16j 




5.49 


± 


0.31 










Lumin. [pb^i] [ 


%] 








These results 




818 


5.27 


± 


0.07 




0.14 


World Average 






5.40 




0.40 







TABLE III: The B^/Be branching ratio for the data based on relative and absolute measurements. 



Method 


B^/Be [%] 


Absolute 


95.98 ± 1.93 ± 1.30 


Cross-feed 


94.64 ± 1.95 ± 1.03 


PDG 2008 


98.36 ±9.16 



Carlo simulations. The cross-feed method is reasonably consistent with the ratio of absolute 
branching fractions. 

IV. PROJECTIVE WEIGHTING TECHNIQUE 

We extract the helicity basis form factors using the projective weighting technique more 
fully described in Ref. |7[. For a given bin, a weight designed to project out a given 
helicity form factor, is assigned to the event depending on its and 6^ decay angles. We 
use 25 joint Acos6'y x A cos 6'^ angular bins: 5 evenly spaced bins in cos^y times 5 bins in 
cos^^*^ 

For each qf bin, we can write the bin populations Ni as a sum of the expected bin 
populations from each, individual form-factor product contribution to Eq. (|2]). Thus 
can be written as a linear combination with coefficients faiQi), 

N^ = h{ql) m+ + rn_ + Uq}) mo 

) rfii + fT{q}) rfiT + hMl) mri • (H) 

Each of the six faiQi) coefficients is associated with one of the form factor products that we 
wish to measure. The six rfia vectors are computed using SPD Monte Carlo simulations gen- 



^ When we use a combined semielectric and semimuonic sample, we use a 50 component N vector with 
the first 25 angular components reserved for £)+ — )■ K~TT^e~^i'e candidates and the second 25 angular 
components reserved for I?+ — > K^-k^ fi^v^ candidates. 
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D-^K+TI-TT- D-^K+ 71-71-71° 




AE (GeV) 

FIG. 5: Distribution of IS.E for single tag D~ candidates when and candidates have been 
combined. The distribution for each of the six tags is shown in (a)-(f). The points with error bars 
are the reconstructed yield from the data sample and the curves show our fit to the signal peak 
over the dashed background line. 

erated with the Eq. ([2]) intensity but including just one of the six form factor products. For 
example, m_|_ is computed using a simulation generated with an arbitrary function for 
(such as Hj^{q^) = 1) and zero for the remaining five form factors. The faiQi) functions are 
proportional to the true H^{qf) along with multiplicative factors such as {Vcsf {.<f — ^j) 
and acceptance corrections. 

Reference shows how Eq. ( ITTll can be solved for the six form factor products H^{q'^), 
H'i{q^), Hl{q^), hQ{q^)HQ{q^), H^{q'^), and HQ{q^)Ht{q^) by making six weighted q^ his- 
tograms. The weights are directly constructed from the six fha vectors. 

Figure [9] shows the six form factor products multiplied by q^ obtained from a Monte 
Carlo simulation using our selection requirements. Because the isolated — )■ K~Ti^e^Ve 
sample provides no useful information on the mass-suppressed form factor products H^{q^) 
and Ho{q'^)Ht{q'^), the second point is not plotted for these two form factor products. The 
Monte Carlo sample was generated with our SPD Monte Carlo with = and was run with 
three times our data sample. The reconstructed form factor products in the Monte Carlo 
simulation are a good match to the input model indicating that the projective weighting 
method is reasonably unbiased. 

We turn next to a discussion of our normalization convention. Equation (j4]) tells us 
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FIG. 6: Distribution of AE for double tag events for the data, where D~ candidate is reconstructed 
in one of the six tag modes [(a)-(f)], and candidate is reconstructed in K*^e'^i^e niode. The 
points with error bars are the reconstructed yield from the data sample and the curves show our 
fit to the signal peak over the dashed background line. 



that as g2 _> o, q'Hl{q') ^ 0; and q^H'^i^l'), q'ho{q')Hoiq'), q'H^{q'), q'Ho{q')H,{q') all 
approch the same constant. Therefore, we normalized the form factor products in Fig. IHlby 
scaling the weighted histograms by a single common factor so that q^Hl{q^) = 1 as — i- 
based on the q'^H^lq'^) measured in the combined — )■ K~7r~^e'^i'e and — )• K~tt^ fi^u^ 
sample. 

Figure [9] shows that the isolated K^ti^ fi^u^ and — ^ K^7c~^e~^i>e samples 

produce similar error bars for the measured H^lq"^), ifi(g^), and ho{q^)HQ{q'^) form factor 
products, while the H^lq"^) errors are much larger for the — )■ K~ t:^ v sample than 
for the — )■ K~'K^e'^Ue. This is due to the large correlation between the Hl[q^) and 
H^{q^) form factors present in the — )■ K^tt~^ fi'^u^ sample owing to the similarity in 
their associated angular distributions. For this reason, the error bars on the H^{q'^) form 
factor product are dramatically reduced when one combines the — )■ K~'k^ ^^v^ and 
— )■ K~TT^e'^Ve samples. 
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AE (GeV) 

FIG. 7: Distribution of AE for double tag events, where D~ candidate is reconstructed in one 
of the six tag modes [(a)-(f)], and candidate is reconstructed in K*^ ^'^u^ mode. The points 
with error bars are the reconstructed yield from the data sample and the curves show our fit to 
the signal peak over the dashed background line. 



V. FORM-FACTOR RESULTS 

We turn next to a discussion of our form-factor measurements. Figure [TD] compares the 
ho{q'^) Ho{q'^) distribution below the nominal pole (a) to that above the nominal K*^ pole 
(b). Figure [TOl shows that there is no significant ho{q^)HQ{q'^) signal above the K*^ pole. The 
absence of a hQ{q'^)Ho{q'^) signal above the nominal K*^ shows that our data are consistent 
with the 6s phase obtained in Refs. 0^ • A related interference pattern was observed in the 
FOCUS [8|] discovery of the s-wave interference in — )■ K~7i~^fi~^h'^ decay. We can thus 
improve our statistical errors by restricting our hQ{q^)HQ[q^) measurements to events with 
0.8 < rriK-K < 0.9 GeV/c^. This additional requirement was applied to the q^hQ{q^)HQ{q^) 
plot of Fig. [m while the other five form factor products use the full 0.8 < niK-n < 1-0 GeV/c^ 
mass range. 

Figure [TT] shows the six form factor products multiplied by q^ obtained for data using our 
q^H^^q^) = 1 as — )■ normalization convention. The background was subtracted using 
our Monte Carlo samples. Although the data are a reasonably good match to the SPD 
model for the g^ifo(g^) and q^H^{q^) form factors, the model does not match the data for 
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FIG. 8: Results on the relative branching ratio, Bfj,/Be obtained for the six tag states and the error 
weighted average of these six values. We compare the relative branching ratio using the cross- feed 
method [Eq. ([6])] to the ratio of absolute branching fractions. Table IIIII gives a summary of these 
results. 



q'^hQ{q'^)Ho{q'^), and the mass-suppressed form factors q^H^{q^) and q"^ Hoi^q"^) Ht{q^) . These 
disagreements will be discussed in Sec. IVTl 

Because of our excellent q^ resolution, there is negligible correlation among the ten q"^ 
bins for a given form factor product, but the relative correlations between different form 
factor products in the same bin can be much larger. Most of the correlations are less than 
30 %. There are, however, some very strong (> 70%) correlations for H_{q^) with various 
other form factors - most notably in the three lowest q^ bins in the correlations between the 
H_{q^) and the Ht Hq as well as Hq form factor products. 

Table IIVI a tabular representation of Fig. [TT] for the — K~tt^ ^'^v^ and — )■ 
K'T^'^e'^Ve combined sample, gives the center of each q^ bin, the measured form-factor 
product, its statistical uncertainty (first error) and its estimated systematic uncertainty 
(second error). The biggest source of the systematic uncertainty is from the background 
estimation. We separately consider systematic uncertainties from non-semileptonic decay 
backgrounds, and semileptonic decay backgrounds. The semileptonic backgrounds include 
cross-feed as well as semimuonic events where the pion and muon are exchanged. 

For the background uncertainty, we assign a conservative systematic error by increasing 
the level of the non-semileptonic background and semileptonic background subtractions by a 
factor of 1.5 and comparing these form factor products to the results with the nominal back- 
ground subtractions. For if^(g^) and hQ{q'^)HQ{q'^), the non-semileptonic and semimuonic 
background subtraction systematic uncertainty is less than 20 % of the statistical error, 
while for the other four form factor products the systematic error is less than 40% of the 
statistical error. 

We also relative systematic error due to uncertainties in track reconstruction 

and particle identification efficiencies. The systematic uncertainty from this source is rather 
small since we are reporting form factor shapes rather than absolute normalization. This 
uncertainty is estimated as less than 1.9 % for all the form factor products. Finally, we assess 
a scale error of 13.4% on the hQ{q^)HQ{q^) form factor product due to the uncertainties in 
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FIG. 9: Non-parametric form factor products obtained for the SPD Monte Carlo sample (multiplied 
by q^) for ten, evenly spaced (f' bins. The reconstructed form factor products are shown as the 
points with error bars, where the error bars represent the statistical uncertainties. The three 
points at each value are: filled circles a combined — )• if ~7r"^//"'"f^j &; — )• K~ tt^ e'^ Ue 
sample, empty squares — t- K^tt^ e'^Ve only, and empty triangles — t- K^tt^ fj.^v'^ only. The 
solid curves represent our SPD model, which was used to generate the Monte Carlo sample. The 
histogram plots are: (a) q^Hl{q'), (b) q^H^{q'), (c) q'Hi(q'), (d) q^ho{q^)Ho{q^), (e) q^Hfiq^), 
and (f) q^Ho{q^)Htiq^). 



the A and 6 values reported in Ref. j9[. When this s-wave scale error is added in quadrature 
to the subtraction systematic error, the total systematic error rises to about 85% of the 
statistical error in the lowest three bins of the hQ{q^)HQ{q^) form factor product, but 
systematic errors on the form factor shape are less than 20% of the statistical error. 

Figure [12] illustrates our sensitivity to the pole masses in Eq. (jl]) by comparing measure- 
ments of the H'^{q^) form factor product to a model with spectroscopic axial and vector 
pole masses versus a model with infinite pole masses, implying constant axial and vector 
form factors. Our data favor the spectroscopic pole masses given in Eq. (|5]), for the high 
bins of the H'^{q^) form factor product. The other five form factor products are consistent 
with either pole mass choice. 

It is of interest to search for the possible existence of additional non-resonant amplitudes 
of higher angular momentum. It is fairly simple to extend Eq. ([2]) to account for potential 

(i-wave or /-wave interference with the K Breit-Wigner amplitude. We search specifically 
for a possible zero helicity d-wave or /-wave piece that interferes with the zero helicity 
K*^ contribution. One expects that such potential h^\q'^) and h\l\q^) form factors would 
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FIG. 10: We show uncorrected plots of the ho{q'^) HqIq^) for data with — K~7r~^ fi'^v^ and 
D'^ — )• K~7r^e~^i'e combined, (a) is for events below the nominal K*^ pole: 0.8 < rriKiT < 
0.9 GeV/c^. (b) is for events above the nominal pole: 0.9 < rriKTr < 1-0 GeV/c^. There is a 
strong ho{q'^)Ho{q'^) signal below the nominal pole but no evidence for a non-zero ho{q^)Ho{q'^) 
form factor above the pole. Note the order of magnitude difference in the y-axis scales between 
the left and right plots. 
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FIG. 11: Non-parametric form factor products obtained for the data (multiplied by q^) for ten 
evenly spaced bins. The reconstructed form factor products are shown as the points with error 
bars, where the error bars represent the statistical uncertainties. The three points at each q^ value 
are: filled circles a combined — )• K~7r~^ n^V/^ & Z)+ — > K^-K^e^Vg sample, empty squares — >■ 
if^TT+e+z^e only) and empty triangles D'^ — > K~'k^ ix^v^ only. The solid curves show our SPD 
model. The histogram plots are: (a) q'^H^q^), (b) q^H'i{q^), (c) q^H^{q^), (d) q'^ho{q'^)Ho{q^), 
(e) q-'HHq''), and (f) q^H^{q^)Ht{q^). 
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TABLE IV: Summary of form factor product results for ten, evenly spaced g^bins for the — ^ 
K^TT^ ^^v^ and — )• K~t:^ e^v^ combined sample. The first and second errors are statistical 
and systematical uncertainties, respectively. The numbers are normalized using the condition: 
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peak as l/y/q^ near — as is the case for the other zero heUcity contributions Ho[q^) 
and holq"^). If so, the zero helicity contributions should be much larger than potential d- 
or /-wave ±1 helicity contributions. The rf-wave projectors are based on an additional 
interference term of the form 

4 sm^ee{3 cos^^y " 1) ^o(g^) hi^\q'^) ReiAae-'^-' (3} . (12) 

To search for the presence of zero helicity d-wave amplitude we use the technique of Ref. 0] 
to construct a projector which is orthogonal to the projectors for each of the six terms in 
Eq. ([2]). The /-wave weights are based on an additional interference term of the form 

4 sin2^^(5 cos^^v - 3 cos ^y) ^o(g^) h[/\q^) Re{Afe-'^f /3}. (13) 

Averaging over the Breit-Wigner intensity, the interference should be proportional to 
A^j sin 6dj h^Q'^\(f) Ho{q'^) and will disappear when the non-resonant amplitude is orthog- 
onal to the average, accepted K*'^ amplitude. Fig. [T^ shows the q^h^Q'^\q^)HQ{q^) form 
factor products obtained in the data using projective weights generated assuming a phase 
of zero. The projective weights are normalized so that q'^h^'^'^\q'^)Ho(q'^) = 1 in the limit 
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FIG. 12: Evidence for finite pole masses. We show the measured q^H^{q^) form factor shown 
in Fig. [11] overlayed with two models, (a) uses the same SPD model shown in Fig. [11] while (b) 
overlays the data with a SPD model where the axial and vector poles [Ma and My in Eq. ^] are 
set to infinity. We show the data with — )■ K^tt^ ^^v^ and — t- K^ir^e'^Ve combined. The 
slight scale difference between the data points in the two plots is an artifact of our q^HQ^q^) = 1 
as g'^ — )■ normalization scheme, which is based on the two different pole mass predictions for the 
HQ{q^) form factor product. 



— )■ if the putative dj -wave amplitude had the same strength as the s-wave amplitude 
relative to the K Breit-Wigner amplitude. There is no evidence for either a (i-wave or 
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FIG. 13: Measurements of the d-wave form factor product (a) and /-wave form factor product 
(b) for an assumed phase of radians relative to the K Breit-Wigner amplitude. 



f-wave component with this phase. 

Figure [H] shows our amplitude and limits for sixteen phase assumptions. As illustrated by 
Fig. [TOl our ability to measure a non-resonant amplitude can depend critically on its phase 
relative to the average, accepted K phase. In order to maximize our sensitivity to the non- 
resonant amplitude, for each phase assumption and bin we made our measurement based 
on three mKir mass regions: 0.8 < < 0.9 GeV/c^, 0.8 < rriKn < 1-0 GeV/c^, and 0.9 < 
mxTT < 1-0 GeV/c^, which puts the average K*'^ reference phase at roughly Sit/A, 3it/2, and 
77r/4 for these three mass regions, respectively. We chose the mass region with the smallest 
expected error according to the Monte Carlo simulation. Under the assumption hQ''^\q^) = 
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FIG. 14: Search for d-wave, (a) and (b), and f-wave, (c) and (d), interference effects for each phase 

assumption as described in the text. The phases 6d and 5f represent the phase of possible d and 

— *o 

f -wave contributions relative to the phase of the K Breit-Wigner amplitude. They are measured 
in radians. 



Hoiq'^), usedinRef. 0, we performed a fit of Fig. [13] to the form Aaj sin 6dj Hl{q^) over 
the region < 0.6 GeV^/c'^ to find the amplitude and limits shown in Fig. UM 

Figure UM shows that this "mass selection" method produced non-amplitude limits, which 
are reasonably independent of assumed phase. If, on the other hand, one used the full 
0.8 < mxTT < 1-0 GeV/c^ mass range for all sixteen phase assumptions, one would get 
dramatically poorer limits for phase choices orthogonal to the Breit-Wigner amplitude phase. 
It is apparent from Fig. [TUthat we have no compelling evidence for either a d-wave, or an 
/-wave component. 



VI. SUMMARY 

We present a branching fraction and form factor analysis of the — )■ K~TT~^i~^i'£ decay 
based on a sample of approximately 11800 — )■ K~7r~^e'^i'e and — )■ K~it^ fi^u^ decays 
collected by the CLEO-c detector running at the ^^(3770). We find BeiD+ K*^e+i/e) = 
(5.52 ± 0.07 ± 0.13)% and Bf,{D+ K*%+u^) = (5.27 ± 0.07 ± 0.14)%. Our direct 
measurement of the relative semimuonic to semielectric branching ratio using Eq. (jG]) is 
B^/Be = (94.64 ± 1.95 ± 1.03)%. 

We also present a non-parametric analysis of the helicity basis form factors that control 
the kinematics of the — )• K~7i~^i~^i'e decays. We used a projective weighting tech- 
nique that allows one to determine the helicity form factor products by weighted histograms 
rather than likelihood based fits. We find consistency with the spectroscopic pole dominance 
model for the dominant if^(g^), ifi(g^) and HQ{q^) form factors. Our measurement on the 
hQ{q'^)HQ{q^) form factor product suggests that the form factor falls faster than Hq with 
increasing q^. The form factors determined using — )■ K~t:^ ix^v^ decays are consistent 
with those determined using — )■ K^Tt'^e^Ve decays and are consistent with our earlier 
study [13] of — )■ K'TT^e^Vf.. Our measured H'^{q^) form factor data are more consis- 
tent with axial and vector form factors with the expected spectroscopic pole dominance q^ 
dependence than with constant axial and vector form factors. 

Our measurements of the H^[q^) and Ho{q'^)Ht{q^) form factor suggests a much smaller 
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Ht{q'^) form factor than expected in Lattice Gauge Theory models Within the context 
of the spectroscopic pole dominance model Eq. ([5]), our HQ{q^)Ht{q'^) measurements are 
most consistent with a small Ht{q^) form factor contribution implying a very negative value 
for r3 = A2,{0) / Ai{0) , such as r^, = —10, which would place the predicted B^jBe relative 
branching ratio close to the phase space estimate of 91%. Finally, we have searched for 
possible d-wave or /-wave non- resonant interference contributions to K^t^^ l^vn. We 

have no statistically significant evidence for d-wave or /-wave interference, but are only able 
to limit these terms to roughly less than 1.0 and 1.5 times the observed s-wave interference 
for d-wave and /-wave respectively. 
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